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ABSTRACT 

The Japan Meteorological Agency started the 
spectral observation of solar ultraviolet (UV) 
ir radiance on 1 January 1990 at Tateno, Aero logi- 
cal Observatory in Tsukuba (36* N, 140* K). The 
observation has been carried out using the Brewer 
spectrophotometer for the wavelengths from 290 to 
325 nm with a 0.5 nm interval every hour from 30 
minutes before sunrise to 30 minutes after sunset 
throughout a year. Because of remarkable simila- 
rity within observed spectra, an observed spect- 
rum can be expressed by a simple combination of a 
reference spectrum and two parameters expressing 
the deformation of the observed spectrum from the 
reference. By use of the relation between one of 
the deformat ion parameters and the total ozone 
simultaneously observed with the Dobson spectro- 
photometer, the possible increase of UV irradian- 
ce due to ozone depletion is estimated. For dama- 
ging UV (I)UV; Wester , 1 98 1 ) f the ir radiance possi- 
bly increases about 19% with the ozone depletion 
of 10% at. noon throughout the year in the nor- 
thern m i d - 1 at i t udos . DU V at noon on the summer 
solstice possibly increases about 5.6% with the 
ozone depletion of 10 m atm-cm for all latitudes 
in the Northern Hemisphere. 

1 . INTRODUCTION 

The possible adverse effects of solar ultra- 
violet. (UV) ir radiance increase due to ozone 
depletion have become of great concern in the 
world. Theoretical studies of UV increase due to 
ozone depletion have been made extensively since 
the 1 970 * s (for example, Cutchis, 1974; Venkates- 
warnn, 1974; Dave and Hal pern, 1976; Pyle and 
Derwent, 1980; Gelstl et al., 1981), whereas the 
studies based on observed UV data have not been 
made enough so far. 

The Japan Meteorological Agency ( JMA ) star- 
ted the spectral observation of global irradiance 
of solar ultraviolet (UV) on 1 January 1990 at 
Tateno, Aerological Observatory in Tsukuba (36.05 
N, 140.00 E). The observations were also started 
at Sapporo, Kagoshima and Naha on 1 January 1991, 
and at Syowa Station, Antarctica, on 1 February 
1991. This paper aims to explain the result of 
the analyses of observations in Tsukuba in 1990, 


focusing on the assessment of the possible UV 
increase due to ozone depletion. The details of 
the entire work have been published in Japanese 
(Ito ct al . , 1991 ) . 

2. OBSERVATION 

The Instruments used in the JMA network are 
the Brewer spectrophotometer Improved on the old 
version used by Josefsson (1986). The important 
improvements are: (1) The polarization prism Is 
substituted by the film polarizer placed on the 
filter wheel, so that the polarizer is only used 
for ozone observations and not for UV spectral 
observation. (2) The width of the spectral scan 
is extended to 290 - 325 nm from the old ver- 
sion’s 290 - 320 nm . (3) The software of the 
Brewer is amended so as to enable the scheduling 
of daily observation based on local time, instead 
of solar zenith angle which is used In the old 
version, (4) The software written with GW- BASIC 
is translated to that with N88-BAS1C to allow 
using a NEC Personal Computer as a control unit, 
which enables us to display the observed spectra 
on the CRT after each observation. 


In Tsukuba, the Brewer #52 is installed near 
the southwest corner of the roof-top of the 
observatory where nothing bars the view field of 
the UV port of the Brewer. The UV port is placed 
0.9 m above the roof floor, 14 m above the ground 
and 39 m above the sea level. The Brewer is ope- 
rated in its automatic continuous operation mode 
and the observations are made every hour from 30 
minutes before sunrise to 30 minutes after sunset 
every day. The measurement of UV spectra takes 
about 8 minutes per once, scanning the wave- 
lengths with a 0.5 nm interval, upward from 290 
to 325 nm, and then downward from 325 to 290 nm. 

The calibration of the Brewer has been con- 
ducted at the Observatory against the NBS stan- 
dard lamp (1000 Watt). From the regular internal 
and external lamp tests, the responsivity change 
less than 7% is noted in 1990. The observed 
spectra arc corrected for this change prior to 
the analyses. 
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3. ANALYSIS OF OBSERVED SPECTRA 

3.1 Parameterization of observed spectra 

The observed spectra show remarkable simila- 
rities in both the general and fine structures of 
the spectra. Noting this similarity, an observed 
spectrum is expressed by a simple combination of 
a reference spectrum at a given solar zenith 
angle and two parameters expressing the deforma- 
tion of the observed spectrum from the reference. 
The reference spectrum is determined as a func- 
tion of solar zenith angle, based on the observed 
spectra in the following procedure. 

Figure 1 shows the scatter diagrams of mea- 
surements with the Brewer (photon count rate) and 
solar zenith angles. The data is taken from those 
for the wavelength channel 305 nm and the period 
from 1 April to 30 June. The curve along the 
upper boundary of plotted data is approximated in 
a form of Eq. 1. 


logCKU, za) =£) OiO-) -cos (it-(i-l) -za/180) (1) _ 

i-1 

where X is wavelength; za is solar zenith angle; 

CR(A , za) Is photon count rate; Q i (XKi = 1.2 8) 

is the coefficient to be determined for each 
wavelength. Note that a scaling of X = (L-325)/ 

0.5 is used to convert the wavelength L in nm 
unit to the nondimensional wavelength 

After the determination of Qj(X) f° r 71 
wavelength channels, the reference UV spectral 
irradiance F(X,za) is determined by 

F a,za)-- CR{l ' za] R ™ {290 ' z ^ (2). 

where R(X) is the spectral responsivity of the 
Brewer #52. The count rate CR(290,za) at the 
shortest wavelength is assumed as the stray light 
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Fig. 1 Scatter diagram of photon count rates 

plotted against solar zenith angles for 
the 305 nm wavelength channels of the 
Brewer, obtained from 1 April to 30 June 
in 1990. 


and used for the stray light correction of the 
spectral irradiance. 

Figure 2 shows the comparison between the 
observed (OBS) and reference (REF) spectrum, and 
the spectral difference (LOR: log( OBS/REL ) ) . Iho 
left ordinate indicates the logarithmic scale of 
spectral irradiance (pW/nT nm), and the right is 
the scale for LOR. 


3.2 Deformation parameter 

The general resemblance of fine structures 
between the OBS and REF spectrum can be seen in 
Fig. 2, which is well demonst rated by the smooth- 
ness of the LOR curve, at least for the wave- 
lengths longer than 300 nm, even though the OBS 
and REF curve exhibits notable fluctuations in 
this wavelength region. The irregular fluctua- 
tions of the LOR curve at the wavelengths shorter 
than 295 nm is due to the low irradiance near the 
detection limit of the Brewer. 

Similar resemblance between OBS and REF is 
seen in almost all observations. Then, the LOR 
curves are approximated for the wavelength longer 
than 300 nm by the following quadratic curve: 


LOR~aX 2 +P 


(3) . 


where a and P are the parameters to be determined 
with each of the observed spectra. 

In Fig. 2. the thick line Indicates the 
quadratic curve thus determined. The fitness of 
the curves can be evaluated by the standard devi- 
ation of the LOR curve around the quadratic 
curve. In all, 4288 spectra are obtained in 1990 
at solar zenith angles ranging from 13 to 90 
degree. Among them, 90% of LOR curves are suc- 
cessfully approximated in the form Eq. 3 with the 
standard deviations less than or equal to 0.06, 
and 1G% are less than or equal to 0.03. The stan- 
dard deviation of 0.03 corresponds to the devia- 
tion of about 1% in irradiance, which is quite 
small compared to the natural variability. 



Fig. 2 Comparison between the observed (OBS) and 
ref erence (REF) spectra of solar U V 
irradiance and the logarithmic difference 
(LOR) between the REF and OBS. 
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Eventually, an observed UV spectral lrradl- 
ance is expressed as 


I(X, za, x , oz, al, •■•) =F(A, za) ‘io aA2+p (4) , 


where t, oz and al are atmospheric optical thick- 
ness (the effect of clouds and aerosols), total 
ozone and ground albedo, respectively. Hereafter, 
a and 0 are called the deformation parameters, 
meaning they express the deformation of an obser- 
ved spectrum from the reference spectrum of a 
relevant solar zenith angle. The parameter a 
Indicates the degree of the wavelength-dependent 
deformation and is expected to correlate mainly 
with oz. The parameter p indicates the degree of 
the wavelength- independent deformation and is 
expected to correlate with r and al . 

3.3 Relation between parameter a and total ozone 

The bottom diagram in Fig. 3 shows the time 
series of total ozone observed with the Dobson 
spectrophotometer for the period from 1 January 
to 31 December 1990. The top diagram in Fig. 3 
shows the time series of the deformation paramet- 
er a at the same time when the total ozone is 
observed, interpolated with the hourly values. 
Comparing these two time series in Fig. 3, it can 
be noticed that the change of a is in good accord 
with that of total ozone. With careful examina- 
tion, however, it is realized that the time 
change of a is perfectly identical In phase to 
that of total ozone, but it is not the case in 


ALPHA (x -10^5) 



TOTAL OZONE Or, atrr-c n) 



amplitude, suggesting that the sensitivity of a 
to total ozone is not uniform, but depends on 
solar zenith angle. 

The total ozone is observed at the solar 
zenith angle between 13' and 70* in Tsukuba. All 
the simultaneous pairs of total ozone oz and 
deformation parameter a are classified into six 
groups based on the solar zenith angle at which 
the total ozone is observed. Each group is chara- 
cterized by the different solar zenith angle 
range with a 10* interval, starting from 10* . 
Figure 4 shows an example of the scatter diagrams 
of oz and a for the solar zenith angles from 60' 
to 70* . The pairs of oz and a in this group 
amounts to 304 and the correlation coefficient 
between oz and a is 99.2%. Similar good correla- 
tion is seen in all solar zenith angle intervals. 

Table 1 shows the mean solar zenith angles, 
the regression coefficient that Is the gradient 
of the regression line; (da/doz)), the offsets of 
the regression lines, standard deviations of a 
(self standard deviations ) , standard deviations 
of a from the regression line (regression stan- 
dard deviations), correlation coefficients, and 
data amounts, corresponding to each solar zenith 
angle interval. 

3.4 Increase of monochromatic irradiance due to 

ozone depletion 

The logarithmic differentiation to Eq. 4 and 
the multiplication of each side by 100 will lead 
to Eq. 5: 


I (a.) doz doz 


(5). 


In this derivation, the dependency of 0 on oz is 
ignored by assuming monochromatic irradiance of 
the maximum wavelength 325 nm is independent on 
the ozone amount. The left side of this equation 


ALPHA ( x-10^5) 



Fig. 3 Comparison of a (top) and total ozone 
with the Dobson (bottom). 


Fig. 4 Scatter diagram of a and total ozone for 
the zenith angles of 60-70' . 
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■eans the percent change of lonochroiatic UV 
Irradiance due to a 1 i ati-cm change of total 
ozone. Hereafter we call it RAC (Radiation Ampli- 
fication Coefficient). 

The values of dct/doz are the regression 
coefficients, listed in Table 1 and plotted In 
Fig. 5 as a function of solar zenith angle. The 
curve in Fig. 5 is expressed by 

da/ doz- const* (cos m {za) +$) “ (1 ** >/in (6) , 


The values of RAC can be calculated with 
Eqs. 5 and 6 for any solar zenith angle, at 
least, between 10’ and 70’ . Figure 6 shows the 
relationship between RAC and wavelength for seve- 
ral solar zenith angles indicated by the numbers 
attached to the curves. It should be noted that 
the values of RAC for the wavelength shorter than 
300 nm In Fig. 6 are given as a extrapolation 
because the values of a are determined based on 
the observed spectra for the wavelength range 
from 300 nm to 325 nm. 


where const * 0.141609, m = 6, 6 = 0.029497 and 6 
= 0.0808, which are determined by the least 
square fitting to the plotted data. 


Tab. 1 Correlation between total ozone and 
deformation parameter a. 

ZA: solar zenith angle; 

GAR: gradient of regression line; 
OFF; offset of regression line; 

ST DEV: standard deviation 
around the average; 

DEV REG: standard deviation 

from the regression line; 
COR COEF: correlation coefficient; 
NUM: number of data. 


ZA 

a 

ST DEV 

DEV REG 

COR COEF 

NUM 

GRA 

OFF 

16.2 

-1.44E-6 

4. 59E-4 

4. 03E-5 

1. 10E-5 

0. 9S18 

28 

25.0 

-1.64 

5.31 

5. 65 

1.44 

0. 9671 

26 

35.2 

-1.68 

5. 74 

5. 95 

1. 17 

0. 9804 

24 

48. 1 

-2.03 

6.71 

6. 87 : 

1.08 

0. 9876 

203 

55.4 

-2.35 

7. 59 

6. 78 

0.92 

0.9308 

79 

66.9 

-2.61 

8.43 

8. 29 

1.05 

0. 9320 

| 308 


d(alpha)/d (oz) 



Fig. 5 Plot of the regression coefficients 

da/doz against the solar zenith angles za 
listed in Table 1. 


3.5 Increase of damaging UV due to ozone 
depletion 

Here, DUV (damaging ultraviolet) is defined 
as the Integration of UV spectral irradiance, 
weighted by the action spectrum formulated by 
Wester ( 1981 ) . The expression of RAC for DUV is 
given by 


loo.dDUV =230 . X 2 da. 
DUV doz doz 


f X2 D(X) * F(X,za ) -10 * x7 **'X 2 dX 



f Xl D(X) 'F(X,za) * 10 * x7 +*dX 


(7), 


( 8 ), 



Fig. 6 Percent increase of the monochromatic IV 
irradiance due to 1 m atm-em decrease of 
total ozone for the indicated solar 
zenith angles indicated by the number 
attached to the curves. 
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where X av is a kind of weighted mean wavelength. 

An assessment of DUV increase due to ozone 
depletion for any geographical location and any 
season can be made using Eqs . 6, 7 and 8, provi- 
ded that the relevant ozone amount and solar 
zenith angle are given. Some results are shown in 
Fig. 7 and in Tables 2 and 3. 

Table 2 shows the monthly DllV-RAC in Tsuku- 
ba. The table also contains the widely used para- 
meter, DUV-RAF (Radiation Amplification Factor) 
indicating the percent change of DUV due to a IX 
change of total ozone. DUV-RAC is computed using 
the monthly average of total ozone and the solar 
zenith angle at noon on the fifteenth day of the 
month. DUV-RAF is obtained by multiplying DUV-RAC 
together with 1% of the monthly average of total 
ozone. The calculated DUV-RAC has the value of 
0.60 ± 0.08 %/m atm-cm , with the minimum in April 
and May, and the maximum in November. DUV-RAF is 
correspondingly given by 1.85 ± 0.2 %/%. Roughly 
speaking, a 10% depletion of total ozone in the 
future will result in a 19% increase of UV ir- 
radiance at noon throughout the year, in middle 
latitudes . 

Figure 7 includes the Northern Hemispheric 
distributions of total ozone, DUV-RAC and DUV- 
RAF. The ozone distribution is drawn on the basis 
of the 11 year averages of TOMS grid data (ver. 
5) corrected by Japanese Dobson values. The value 
of DUV-RAC is computed each grid with the total 
ozone at the grid point and the corresponding 
solar zenith angle at noon of the summer solstice 
(the values in the region between 25* N and 35* N 
are extrapolated). The value of DUV-RAF is com- 
puted each grid by multiplying RAC by 1% of the 
total ozone. In Fig. 7, It can be seen that the 
longitudinal variation in RAC (and also RAF) is 
less remarkable compared to the latitudinal vari- 
ation. 

Table 3 shows the longitudinal averages of 
the grid data shown in Fig. 7. The value of DUV- 
RAC marks the minimum between 45* N and 50* N, and 
the maximum at the equator and at latitudes bet- 
ween 75* N and 80* N. On the average, DUV-RAC 
ranges within 0.555 ± 0.035 %/m atm-cm, and DUV- 
RAF within 1.81 ± 0.28 %/%. It can be concluded 


Tab. 2 Radiation amplification coefficients 

(RAC) and radiation amplification factors 
(RAF) of damaging ultraviolet (DUV) for 
noon in the northern mid- iat I tudes each 
month (MON). 


MON 

Z A 

O Z 

RAC 

RAF 

1 

57.2 

331 

0 62 

2.05 

2 

48.8 

351 

0 58 

2. 05 

3 

38. 2 

362 

0 54 

1. 95 

4 

26 3 

354 

0.52 

1 . 86 

5 

17.2 

349 

0. 52 

1.82 

6 

12.8 

337 

o r.3 

1.78 

7 

14.5 

313 ! 

0 54 

1.70 

8 

22.0 

294 

ij. 50 

1 65 

9 

33 0 

287 

0. 59 

1 70 

10 

44. 5 

280 , 

0 64 

1 80 

11 

54.5 

282 ; 

0. OH 

1 .92 

12 

59.3 

307 

0. 65 

1 . 98 


OZONE 



RADIATION AMPLIFICATION COEFFICIENT* 1000 



RADIATION AMPLIFICATION FACTOR*100 



Fig. 7 Northern hemispheric distribution of 
total ozone (top), the radiation 
amplification coefficient (middle; 
multiplied by 1000) and the radiation 
amplification factor (bottom; multiplied 
by 100) at noon in mid-summer. 


Tab. 3 Radiation amplification coefficients 

(RAC) and radiation amplification factors 
(RAF) of damaging ultraviolet (DUV) at 
noon in mid-summer of each latitude 
(LAT) of the Northern Hemisphere. 


LAT 

Z A 

' OZ 

RAC 

RAF 

0 

23 3 

i 261 

0. 59 

: i 53 

5 

18.3 

268 

0. 58 

i 1.54 

10 

13.3 

275 

0 57 

1 1 57 

15 

8.3 

281 

0.57 

1 1.60 

20 

3.3 

288 

0 57 

1.63 

25 

1. 7 

295 

0. 56 

I. 66 

30 

6.7 

302 j 

0. 56 

1. 67 

35 

11. 7 

313 : 

0. 55 

1. 70 

40 

16. 7 I 

331 

(! 53 

1 76 

45 

21. 7 

351 i 

0 62 

; 1 63 

50 

26.7 

361 

fj 62 

1 1. 88 

55 

31. 7 

363 1 

0. 63 

! 1 91 

60 

36.7 

350 

i; 54 

1 . 94 

65 

41. 7 

353 

n 56 

1 97 

70 ! 

46. 7 i 

353 

11 57 

2 03 

75 

517 i 

356 

0 f.'i 

2 OH 

80 

56. 7 1 

353 • 

0 59 

3 09 

8S 

61. 7 

355 | 

(]. 67 

2 02 


661 



from Table 3 that DUV at noon on the summer sol- 
stice would increase by about 5.6% due to the 
ozone depletion of 10 m atm-cm throughout the 
Northern Hemisphere. 

The present method to evaluate RAF or RAC 
for DUV can be similarly utilized for other ac- 
tion spectra. It may also be used for assessing 
the geographical distribution of UV increase due 
to the predicted geographical distribution of 
ozone depletion. 

4 . SUMMARY 

The procedure estimating the UV increase due 
to ozone depletion is described, based on the 
simultaneous observation of spectral global UV 
irradiance with the Brewer and the total ozone 
with the Dobson, Noting the similarity in shapes 
among observed spectra, the reference UV spectrum 
is determined as a function of solar zenith 
angle, based on the observed data, which give a 
practical maximum of spectral irradiance for a 
given solar zenith angle. An observed spectrum is 
expressed by a simple combination of the referen- 
ce spectrum and two parameters expressing the 
deformation of the observed spectrum from the 
reference. The percent Increase of monochromatic 
UV irradiance due to the ozone depletion of 1 m 
atm-cm is evaluated by using the regression coef- 
ficient between the deformation parameter and 
total ozone, determined as a function of solar 
zenith angle. Also, the seasonal and geographical 
distribution of the percent change of DUV due to 
ozone depletion is evaluated. 

The results show the DUV irradiance inc- 
reases about 19% with the ozone depletion of 10% 
at noon throughout the year in the northern mid- 
latitudes. Furthermore, the DUV at noon on the 
summer solstice increases about 5.6% with the 
ozone depletion of 10 m atm-cm for all latitudes 
in the Northern Hemisphere. 
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